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Abstract. In 2002 a measurement of the effect of solar gravity upon the phase of 
coherent microwave beams passing near the Sun has been carried out with the Cassini 
mission, allowing a very accurate measurement of the PPN parameter 7. The data have 
been analyzed with NASA's Orbit Determination Program (ODP) in the Barycentric 
Celestial Reference System, in which the Sun moves around the centre of mass of 
the solar system with a velocity Vq of about 10 m/sec; the question arises, what 
correction this implies for the predicted phase shift. After a review of the way the ODP 
works, we set the problem in the framework of Lorentz (and Galilean) transformations 
and evaluate the correction; it is several orders of magnitude below our experimental 
accuracy. We also discuss a recent paper |10j . which claims wrong and much larger 
corrections, and clarify the reasons for the discrepancy. 
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1. Introduction 

In 2002 the Cassini spacecraft, on its cruise to Saturn, has allowed an outstanding 
measurement of the effect of solar gravity on the phase of a coherent microwave beam 
sent from the ground antenna to the on-board transponder and transmitted back. The 
data have been analysed with the Orbit Determination Program (ODP), a numerical 
code developed at NASA's Jet Propulsion Laboratory for accurate space navigation, 
including relativistic effects; its algorithms are described in detail in [H] (referred to 
here as the Manual). In the rest frame of the Sun the total phase change can be measured 
by the light-time t2 — ti between the events 1 and 2, the start and the arrival of a photon 
in the up- or the down-link, respectively: 



Here ri and r2 are the distances of 1 and 2 from the Sun at their appropriate times ti and 
h', fi2 is their Euclidian distance; m = 1.43 km is the gravitational radius of the SuiJjI. 
This formula shows a characteristic enhancement in the delay At near conjunction, 
when the impact parameter b of the ray is much smaller than both ri and r2] in this 
case the argument of the logarithm in ([1]) is about 4rir2/6^, as in eq. (1) of [4J. This 
approximation, however, is not accurate enough and is never used in the ODP. The 
fractional change in frequency Ay is essentially the rate of change of the light-time; as 
one can intuitively see, and as fully discussed in [3j (quoted as BG)|§|, this frequency 
shift is of order av, where a is the deflection angle and v a linear combination of the 
velocities of the end points. The fit of the Cassini data gave the result: 



An important feature of the experiment was the use of a multi-frequency link, which 
allowed an excellent elimination of the contribution of the plasma in the solar corona. 
To date, this is by far the best measurement of the PPN parameter 7, which has a 
crucial role in discriminating between alternative theories of gravity. 

Due to the other planets, in particular Jupiter, the Sun moves around the centre 
of gravity of the solar system with a velocity f ~ 5 x 10~^ (15 m/sec) and a time scale 
of several years. During a light-time (about an hour) its motion is essentially uniform, 
so that its effect on the light-time can be described by a Lorentz transformation and 
has nothing to do with the intrinsic scattering dynamics. The short Nature paper [1] 
does not mention the motion of the Sun. The Manual does not address this problem 
explicitly, nor evaluate the correction; however, we show that the ODP does take into 
account this effect, and the fractional correction to the one-way gravitational delay At 

X We mainly follow the notation of [4 and the Manual [M]; the velocity of light, however, is c = 1. 
Eq. (1) of [4j refers to the round-trip and has an additional factor 2; in eq. (2) an obvious factor 2 is 



§ G. Giampieri died on Sept. 2, 2006 at the age of 42. We take this opportunity to acknowledge his 
outstanding and varied scientific contributions, in particular to BG. 



t2-ti= ri2 + At = ri2 + (7 + l)mln 



ri + r2 + ri2 
ri + r2- ri2 



(1) 




(2) 



missing. 
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is of order Vq, quite below the experimental sensitivity. The effect can also be described 
in terms of the induced fractional frequency shift Ay = av, which attained the value 
3 X 10~^° in the best passage. We show that the motion of the Sun causes a change 
in Ay of order avvQ, much smaller than the frequency noise cXy ~ 10~^^. In a two-way 
experiment, like Cassini's, even this small correction largely cancels out. In Sec. 4 we 
discuss the recent paper [10], which claims wrong and much larger corrections, and point 
out the reasons for the discrepancy. Einstein's prediction 7 = 1 is still unchallenged. 

2. Using the Orbit Determination Program 

The theoretical discussion in the main text of [1] (eqs. (1) and (2)) was given only 
to help physical understanding. The Block 5 receivers at the ground station count the 
cycles of the phase of the received coherent beam relative to the local frequency standard 
and determine the number of cycles in successive time intervals. The duration of these 
intervals can be chosen from 0.1 sec to several hours (Sec. 13.3.1.3 of the Manual); during 
Cassini's experiment it was 1 sec, but further averaging was done in the analysis. The 
extraction of a time derivative from the actual phase count, which contains fast- varying 
components, would be very delicate; the frequency is never measured, nor is needed, 
since the mathematical expression of the predicted light-time is at hand. This is why 
the ODP deals only with light-time. The difference between Cassini's experiment and 
measurements of the radar delay does not lie in a different gravitational observable, but 
in the fact that Cassini uses a coherent radio beam and its phase as the main observable, 
while planetary radar determines the arrival time of a wave packet through the peak of 
its intensity. The strongest useful signal in Cassini's experiment occurred on day 2 after 
conjunction, when the logarithm in ([1]) was about 10, with a delay 

At = 1.43 (7 + 1)10^ cm; 

thus the formal accuracy 2.3 x 10^^ in 7 corresponds to an accuracy of 30 cm in At. 

The ODP includes a very extensive and well tested orbit determination program, 
with which the orbits ri{t) of the centres of all relevant bodies in the solar 
system, including the Sun and any spacecraft, are determined numerically from 
previous and current observations. Following standard astronomical usage and lAU 
recommendations, in particular the 2000 lAU resolutions, this is carried out in the 
solar system Barycentric Celestial Reference System (BCRS). The excellent paper [22] 
(referred here as S) describes the procedure in detail. The expression S(8) for the metric 
tensor is the basis for the definition of this frame; but, as T. Damour has discussed in 
the first [6] of four papers on relativistic frames, strictly speaking the BCRS cannot be 
described in a single coordinate chart, but must use a global chart for the dynamics of 
the gravitating bodies, and a local chart for each gravitating body. In the linearized 
approximation we use, however, this is not needed and S(8) has a global meaning. In the 
BCRS the origin is at the centre of gravity C of the solar system; its axes are anchored 
to the non rotating 'celestial sphere' - the International Celestial Reference System 
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(ICRS) - as realized by very distant astronomical sources. Harmonic coordinates are 
used: to the required order, gravity is described by a ('small') spacetime tensor h^^, 
in a fixed and conventional Minkowsky background with metric 7]^^^ = diag(— 1, 1, 1, 1). 
In this view one can use Lorentz formalism; the word 'null', the proper time r and 
raising and lowering spacetime indices refer to this metric. In addition, the (cartesian) 
space components hmn are proportional to the unit matrix 6mn (the coordinates are 
'isotropic'). The time coordinate is the Barycentric Coordinate Time (TCB), related to 
the Earth-bound Geocentric Coordinate Time (TCG) by eq. S(58), in which, essentially, 
the Doppler shifts due to the motion of the Earth and the gravitational shift due to 
external gravitating bodies are corrected for. 

In the present paper we only need to consider, besides the BCRS, the heliocentric 
reference frame, neglecting quadratic terms in its gravitational potential and the 
contributions of other masses. Then in isotropic coordinates both hoo and the space 
components /imn are functions only of the distance R from the origin. Note that in 
general relativity the field equations do not uniquely determine the metric tensor in 
terms of the radial coordinate; in our special coordinates this gauge freedom is forfeited. 
The time coordinate T is the privileged and invariant Killing time, with respect to which 
the metric tensor is constant. The standard textbook expression for the light-time (see, 
e.g., [25]) reads 

T,-T, = R,, + AT = Ru + (7 + l)mln ^ (3) 

The coordinate^ of the end events are Xi = (Ti, Ri), X2 = (T2, R2). It should be noted 
that this formula holds only in isotropic coordinates. 

The ODP instead (see Fig. 1) computes the light-time in the BCRS, where the 
centre of the Sun rQ(t) moves around C with a velocity f ~ 3 x 10^^; squares of Vq can 
be neglected. As described in Sec. 8.3.6 of the Manual, a special subroutine implements 
the formula in 32 steps, an inclusion discussed at length in [8]. The main object 
of the present paper is to evaluate the difference in the light-time in the two frames. 

In the Manual, Sec. 8.3.1.1, the three quantities ri, r2 and ri2 are defined in 
barycentric coordinates as follows. The numerical code provides for the bodies 1 and 2 
and for the Sun the orbits ri{t), r^(t), TQ{t) relative to C in the BCRS. The simultaneous 
differences (eq. (8-62) of the Manual) 

ri(ti) = rf (ti) - rg(ti), r2(t2) = r^ih) - ^it^) (4) 

provide ri = |ri(ti)| and r2 = |r2(t2)|; is the modulus of 

ri2 = r2(t2)-ri(ti). (5) 

II We use small letters with greek indexes from to 3 to denote spacetime quantities as geometrical 
objects in a generic coordinate system; when time and space (boldface) components are written 
separately, we use special coordinates: small latin letters in the barycentric frame and capital letters 
in the rest frame of the Sun. 

*li This corresponds to eq. (8-55) of the Manual, which is the one actually coded in the program. 
Besides the Sun, (8-55) may include the gravitational delays due to 10 other bodies. The argument of 
the logarithm in ([1]) includes a small correction of order m which is not relevant to our problem. 




Figure 1. The geometry of the experiment in the BCRS frame at epoch June 23-th 
2002 20:00:00 UT, the first useful passage after conjunction in the Cassini experiment. 
The center of the Sun moves in the ecUptic plane around the barycenter of the solar 
system (SSB) in a roughly circular orbit with radius 7.8 x lO^m; its velocity Vq = 14.6 
m/sec is approximately orthogonal to the direction of Jupiter, and the orbital period 
is « 12 y 

Note that differences between positions of bodies taken at the same time are invariant 
under Galilei transformations; hence ri, r2 and ri2 have, up to 0{vq), the same value 
in the rest frame and the barycentric frame. We'll return to this important point later. 

It is useful to clarify the logical status of eq. ([1]) and of the subroutine used in the 
ODP to implement it, as explained in Sec. 8.3.2 of the Manual. Given the arrival time 
^2 at the point r2(t2)5 the light-time solution must provide the starting time ti, and the 
position and velocity of the transmitter 1 at that time. A first approximation to ti is 
given by the geometrical solution 

ti=t2- ruiti), 

where the argument of the distance indicates the time at which 1 is taken; an error 
At ^ 10^^ sec is made due to the neglected gravitational delay. In this short time 
interval the velocity of 1 is uniform and a linear correction 6ti 
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ti + St^ = t2- ru{t^) - 5t^-^ - At 

is sufficient. The corrected value of the hght-time is then compared with the 
observations. 



3. Measuring the light-time in the barycentric frame 

We ffist comment about coordinates and geometrical objects. In the ODP scheme - the 
radial coordinate r being fixed - physical laws and statements are invariant under the 
Lorentz group; to deal with this requirement we take the geometric point of view, in 
which scalars, vectors and tensors - like h^^, - denote a geometrical object, independent 
of its coordinate representation. Measured quantities are invariant scalars. The choice 
of coordinates is conventional and free; for a vector, like = {vq, v), they are indicated 
in round brackets, with the space components in bold type. 

We only need to deal with the spacetime vector from 1 to 2, the geometric object 
which summarizes the scattering dynamics of the experiment. In the rest frame, where 
the Sun sits at the origin, its components are 

r = - = (T2 - Ti, R12) = (i?i2 + AT, R2 - Ri); (6) 

if the mass m vanishes 



r = (/2l2,Rl2) 

is null; gravity just adds a small, time-like and positive contribution (AT, 0) along the 

time coordinate. An observer with (constant) spacetime velocity vector measures 

the hght-time —v^i^; if the observer is at rest relative to the Sun, Vq = (1, 0) and the 

previous result T2 — Ti is recovered; if it is at rest relative to the barycenter C, the ODP 

value t2 — ti is obtained. Of course, the components of the light-time in the two systems 

are related by a Lorentz transformation. 

[To] . [9] and other papers take a different approach. They apply the appropriate 

Lorentz transformation to the full trajectory of the photon in the rest frame (as given, for 

instance, in [25]) and obtain its expression in the barycentric frame; then they compute 

the elapsed time. Delay and bending of the transformed orbit are then determined by the 

gravity of a moving Sun, with new gravitomagnetic terms appearing. This complicates 

matters. A similar situation would arise in the classical Rutherford scattering of an 

electron by a proton, normally described by the Coulomb field of the latter, assumed 

at rest; of course, one would get the same result in a frame where the proton moves, 

provided the magnetic field so generated is taken into account. Just as magnetism is the 

direct consequence of the fact that under Lorentz transformations the electromagnetic 

potential behaves as a four-vector, so gravitomagnetic metricperturbations arise due to 

Lorentz invariance and the tensorial character of the metrictU. 

For an elementary derivation of gravitomagnetism from Schwarzschild linearized solution, see [2], P- 
571. 
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The analysis of Cassini's experiment has been carried out in the standard framework 
of a Lorentz invariant theory. If Lorentz invariance is violated, gravitomagnetism does 
not have the standard form; moreover, the problem of the effect of the motion of the 
Sun on the gravitational delay must be addressed in a way different from the one we 
follow below. KPSV suggest that, with its excellent accuracy, Cassini's experiment 
may set better limits to such violations. This would require a wise decision about 
the best theoretical scheme and the appropriate parametrization; besides the PPN 
preferred frame formalism, one has vector-tensor theories (see |25]) and Lorentz- violating 
electromagnetism, on which there is a wide literature (see [15], [16], [18], [Tj). Current 
limits from other experiments should also be taken into account. This program is outside 
the present paper. 

It is more convenient to skip the formal coordinate change in the trajectory and to 
use just the light-time vector - a geometrical object - expressed as 

r = (t2-tl,X2-Xi) 

in the barycentric frame and as 

r = (T2-Ti,R2-Ri) 

in the rest frame of the Sun. The barycentric components are related to ([6]) by a Lorentz 
transformation corresponding to the (small) velocity of the Sun v©: 

t2-ti =r2-ri + Vo-Ri2 + 0(t;^), (7) 

X2-X1 =Ri2+Vq(T2-Ti)+0(4). (8) 

The space part is just a Galilei transformation, arising from the fact that the end 
events occur at different times. In Galileian relativity time is invariant; but this holds 
only if the space component of a vector upon which the transformation operates is 
much smaller than the time component. In our case £^ is almost null and ([7j) must 
be taken into account. The transformation conserves proper length to the appropriate 
approximation: 

ih - t,f - |X2 - xip = (T2 - T,f - Rl, + 0(t;|); (9) 
or, neglecting squares of the delay. 

At |x2 -xil = ATi?i2. (10) 

Now 

X2-X1 = r^(t2)-rf(ti), 

the space component of in barycentric coordinates, obviously is not Galilei invariant; 
but the Manual shrewdly uses ri2 instead, which, being built with simultaneous vectors, 
almost compensates the Galilei transformation ([8]). Indeed, from ([5]) 

ri2 = X2 - xi - V0(t2 - ti); (11) 
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comparing it with ([H]), we see that 

ri2-Ri2 = 0(t;|). 

Let n = ri2/ri2 be the (Gahlei invariant!) unit vector along the unpertubed ray from 1 
to 2, as computed in the Manual. Eq. (IHl) gives, to 0{vq), 



|x2 - xi| = i?i2 + V0 • n(T2 - Ti), 

so that 

At = AT(1 - Vq ■ n), (12) 

which is our final result. It corresponds to a negligible change of about a millimiter 
in light-time. On p. 8-28 the Manual acknowledges the neglect of this correction: The 
error in the calculated delay due to ignoring the barycentric velocity % of the gravitating 
body has an order of magnitude equal to the calculated delay multiplied by the velocity 
of the body/c. Our procedure clearly shows the gist of the problem, how the delay 
transforms under a Lorentz transformation with a small velocity. As C. M. Will has 
pointed out to us, a simpler proof is based upon his eq. (7.3) in [25], which says how the 
perturbation in the coordinate of the photon along the ray changes with time; in the 
barycentric frame (7.3) acquires an additional term due to the time-space components 
of the metric and proportional to m/r; this leads to our result. Note also that in a 
two-way experiment, like Cassini's, in the down-link n is almost equal and opposite its 
up-link value and the two contributions essentially cancel out; a net result arises from 
second order terms 0{Vq) in the Lorentz transformation. 

Eq. ( TT2i) does not say, however, how to evaluate the delay in terms of the distances 
ri,r2 and ri2 obtained in the BCRS. We need their expressions in a generic frame, 
such that they reduce to Ri, R2 and R12 in the rest frame. In the present context, 
invariance under the approximate Lorentz group ([7]) and ([8]) is sufficient. Then the ODP 
expressions (jlj) and ([5]) provide the straightforward answer: in fact they are constructed 
with simultaneous differences and ([H]) reduces to Galilei transformations. 

It interesting to discuss the fully invariant case. Distances arise here because the 
Sun affects a photon essentially through the scalar potential U, the fundamental solution 
of D'Alembert's equation. In the rest frame and, say, for a photon at the start event 
= (Ti,Ri), 

|rv,i| iti 

We want to express it in a generic frame, with the Sun moving with an arbitrary (but 
constant) four- velocity Vq = dX'^/dr. The cartesian distance Ri must be replaced 
with the unique quantity which (i) is Lorentz invariant and (ii) reduces to Ri when 
Vq^ = (—1, 0) (the minus sign being required by the metric diag(— 1, 1, 1, 1)). This is a 
textbook problem in electrodynamics (e. g., see [H], Sec. 63). Let be the (ante- 
dated) intersection of the (straight) world-line of the Sun with the past light-cone of 
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the start event = + £i ; ii is the corresponding, future-pointing null vector. The 
required 'distance' is the invariant 

= \ve,{x't - = Ml . (14) 

Indeed, since 



X 



?-X? = |xi-Xi| 



(elapsed time is equal to distance), it obviously fulfils (ii). The future light-cone would 
give the same result. 

In the appropriate slow motion approximation Vq^ = (— 1, v©) and (fT4l) reads 



X 1 



X°-Vo-(xi-Xi) + O(t;0)2. (15) 



Xi — X°, the antedated distance between the event 1 and the Sun, differs from the 
simultaneous value by an amount equal to the distance traveled by the Sun in a time 
equal to the distance itself; hence is, to 0{vq), just the simultaneous distance ri 
(see dl]) used in the ODP. In other words, the first-order correction to Lienard-Wiechert 
potential due to a slowly moving source vanished. 
The fully invariant form of ri2 = |r2 — ri| is 

ru = \ve,^'\ = \veM - (16) 

indeed, in the rest frame Xi = (Ti, Ri) and Xi = (T2, R2), with T2 — Ti = IR2 — Ri| = 
i?i2. We now show that in the slow motion approximation 

ru = - X? - Vq ■ (X2 - xi) + Oivef (17) 

reduces to modulus of the ODP ([5]). We must use the simultaneous relative vector 
ri = Xi — X(x'j'), which differs from the space part of ii 

^1 = xi - Xi = xi - X(x'|' - ri) = ri + riv© 

because of the ante-dated time argument in the position X(a;°) of the Sun. With a 
similar first-order expansion we get 



X2 - Xi = X(x^ - r2) - X(x? - ri) = (x^ - x[ - r2 + ri)vQ = {r^ - r2 + ri)v0, 
where we have used 

X(x°)-X(x?) = Vo(x°-x?), 

X2 — Xi = |x2 — Xi I and the fact that X2 — Xi differs from r2 — ri = ri2 by 0(fo)^- Then 

* For the full expansion of the retarded potential in powers of 1/c, in which the term 0(l/c) is missing, 
see [7], supplementary note 10 in the Appendix. 
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|x2 - xi| - ri2 ■ Vq = |ri2 + X2 - Xi + (r2 - ri)vQ| - ■ 
|ri2 + ri2V0| - ri2 ■ Vq = + 0{vqY. 



We see that the ODP variable ri2 is the appropriate approximation of a fully invariant 
quantity and describes correctly the gravitational interaction. 

4. About Galilei invariance 

KPSV make two different claims about the effect of the motion of the Sun. Their 
expression of the gravitational delay in the barycentric frame KPSV(25) 

A J = (7 + l)m In ri + r2 + r-r..^ _ 
ri + r2 - r + r ■ V0 

differs from our analysis in two ways. First, it is based upon the coordinate 
transformation KPSV(ll), which does not provide the required, first order difference 
between t and T (eq. ([7])); this makes it impossible to get the true correction f|T2l) . More 
important, the quantity ri2 used in the ODP is not their r, which, in our notation, is 
|x2 — Xi|. This quantity is not a Galilei invariant and cannot appear in the argument 
of the logarithm. This argument must depend on the coordinates of the photon relative 
to the Sun. Let us see the consequences of this claim. In the extreme conjunction 
case {b ^ (^"1,^2)) the correction in the denominator of the argument of the logarithm 
prevails and ( |T8l) reads 



A^t = (7 + l)m [\n - V0 • n-^ j . 



The claimed correction ^ ttivq is similar to the true correction f|T2l) . but is enhanced 
when b is small. 

In their second claim, KPSV consider the description of the experiment in terms of 
the effect that the deflection has on the frequencies recorded by two very distant clocks 
(see [3|). Of course the deflection angle 

Tfl 

a = 2(7 + l)y (19) 

is meaningful only in the extreme conjunction approximation b <^ (^"1,^2), which we 
assume here for the sake of illustration. This gives also, in order of magnitude, the 
change in the angle between the velocity V of a far-away clock relative to the Sun and 
the beam, and corresponds to a (one-way) fractional frequency change (for a grazing 
ray and the Earth's orbital velocity) of order 

Ay^aV^Sx 10^^°. (20) 

Explicitly (see eq. (22) of BG), in the rest frame we have 

^ (V2 ■ b)ri - (Vi ■ b)r2 2(7 + l)m ^ ^^i) 
ri2 b 
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In the following, for simplicity, we only consider frequency measurements in terms of 
the proper times of 1 and 2. Now Ay, a ratio of proper frequencies, is a Lorentz- 
invariant scalar; to get its expression in a slowly moving frame, note first that the time 
coordinate is not involved, so that only Galilei transformations need to be considered. 
We just need the unique quantity which (i) is Galilei invariant and (ii) reduces to (l2Til 
in the rest frame. The impact parameter vector is the component of r2 (or, equivalently, 
— ri) along the ray: 

b = r2-(n-r2)n. (22) 

This is invariant. Secondly, and trivially, the rest-frame vectors Vi and V2 should be 
replaced with their invariant forms 

dri{ti) dv2{t2) , . 

dt ' V2 - V0 = . (23) 

KPSV claim that in the barycentric frame, due to the motion of the Sun, there is 
an additional correction 

Axy ^avQ^Sx 10"^^ (24) 

The fact that the impact parameter vector (parallel to the deflection vector a) is a 
Galilei invariant seems to be ignored; it is defined in the rest frame in terms of absolute 
coordinates, with no reference to the Sun. Then in the barycentric frame the correction 
KPSV(18) arises, which produces a small correction 0{vq) in the frequency shift. 
Moreover, absolute velocities are used in the analogue of f l2T]) . which, after a Galilei 
transformation, produces KPSV(21), the basis for their main claim. Of course, the 
claimed correction (!24l) is absorbed in (12T|) using the relative velocities vi — v©, V2 — v©. 

In Cassini's experiment the standard deviation of the residual frequency noise (at 
an integration time of 1000 sec; see Fig. S2 of the supplementary material^ of [1]) 
was ay = 0.7 x lO"^''. While the magnitude of the correction IHM stated by KPSV 
is 40 times larger, the paper does not even suggest a violation of general relativity. 
In the experiment = 1094 data points, obtained after several pre-processing steps, 
have been used, with a remarkably Gaussian distribution of the frequency residuals, as 
impressively shown by Fig. 3 of [4j. The sensitivity of 7 from each data point is difficult 
to ascertain, but surely the final error in 7 which results from the fit is smaller than 
ay. KPSV, instead, state the opposite (top of p. 280) and thereby avoid dangerous 
conclusions. In our view, if the ODP dynamical model were as inaccurate as claimed by 
KPSV, the inescapable conclusion would be that the good orbital fit and the agreement 
with the standard prediction are the result of an exceedingly unlikely chance; accepting 
the KPSV correction would imply the first experimental violation of general relativity. 

tt http:/ /www. nature. com /nature/journal/v425/n6956/suppinfo/nature01997.html. 
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5. Conclusion 

Cassini's experiment, by far tlie most accurate attempt to measure tlie PPN parameters, 
did not pose a threat to Einstein's theory of gravitation, nor answer the major question, 
at what level and how it is violated. The measurement of the deflection parameter 7 has 
a crucial importance; although no theoretical prediction is at hand about its deviation 
from unity, alternative theories based upon a long range scalar field require 7 < 1. In the 
future steady improvements in instrumentation, in particular space astrometry and the 
use of optical links in interplanetary space, will allow much more accurate measurements, 
and new relativistic space missions are under construction or in the planning phase; see 
[23] and [21] for reviews. We quote, in particular, the astrometric mission GAIA [12], 
the experimental program MORE [13] in the BepiColombo mission to Mercury and 
the project LATOR [23]. Optical interferometric measurements, similar to those under 
development for gravitational wave detection (the LISA project), are planned for the 
ASTROD mission ([IT], [19], [20]). Binary pulsars, whose gravitational field is much 
stronger, will also play a role in testing gravitational theories. 

A numerical code for relativistic orbit determination more accurate than the ODP 
is surely needed. Aside from the correction ( fT2l) in the light-time, terms in the 
metric quadratic in the gravitational radius m must be included. The second order 
gravitational deflection depends on the radial gauge (see |5j); in the isotropic gauge, for 
a grazing ray, it has the value Oquad = 1.2 x 10~^^, corresponding to a frequency shift 
~ foaquad = 10~^^, an accuracy easily accessible with laboratory standards. We should 
also mention that in conjunction experiments, besides the obvious smallness parameter 
m/i?0 = 2 X 10~^, the smallness of b/r may produce an enhancement in the light-time. 
At the first order in m this is already apparent in the logarithmic term of ([T]); but there 
is also a quadratic correction to the light-time of order which may produce an 

effect in Cassini's experiment about the same as the formal error. A correction of this 
kind is included in the ODP. Finally, the complex procedures used in the definition and 
the construction of the BCRS and the TCB need to be revisited and improved both at 
the fundamental level and in the numerical implementation. 
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